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Abstract: The correlation dimension and /^-entropy are estimated from meteorological time-series. 
The results lead us to claim that seasonal variability of weather is under influence of low dimensional 
dynamics, whereas changes of weather from day to day are governed by high dimensional system(s). 
Error-doubling time of this system is less than 3 days. We suggest that the outstanding feature of the 
weather dynamics is deterministic chaos. 

Keywords: correlation dimension - ^-entropy - weather predictability - chaotic dynamics 

1. Introduction 

Prediction and predictability are matters of fundamental importance to meteorology. Atmosphere is a 
considerably complex system containing a wide range of processes linked together by complicated non-linear 
feedbacks. Some people regard weather forecasting as a main task of meteorology (atmospheric physics). One 
can watch on TV weather "forecasts" one month or even several seasons ahead. Hardly anybody of laymen asks 
about success of such long-range prediction. Weather forecasters know that there are meteorological regimes that 
are more predictable than the others. It was already known to John von Neumann in fifties. In conference in 
Princeton on the Dynamics of Climate (Pfeffer 1960) he divided the motions of the atmosphere into the three 
groups: 

1. Processes which are determined by the initial conditions - we may extrapolate the initial trends 
over a short period. 

2. Processes which are practically independent to the initial conditions - if we want to forecast such 
processes we concern with particularities of the atmospheric circulation which will always be 
present. 

3. Between the two extreme cases there is another category of processes - if we are sufficiently far 
from the initial state so that the details of the initial conditions do not express themselves clearly in 
what has developed. Therefore we can't calculate backwards the initial conditions from the ultimate 
conditions. 

In our paper, we focus on using time-series of values of air temperature and atmospheric pressure to determine 
whether or not the atmosphere is chaotic and extent to which atmospheric motion is predictable. The data include 
the daily mean surface temperature and surface pressure of air over period of 30 years (1 January 1961 - 31 
December 1990) observed in the eight Czech meteorological stations (Cheb, Praha - Ruzyne, Hradec Kralove, 
Kuchafovice, Praded, Holesov, Mosnov, Lysa hora). An approach based on technique developed from 
deterministic chaos theory is discussed. 

2. State Space Reconstruction 

The key tool used for identification of chaos in time-series is embedding scheme (Packard et al. 1980, 
Takens 1980, Sauer et al. 1991). This procedure enable us to reconstruct the essential features of dynamics in m 
dimensional state space from scalar time-series. Let us assume that we have recorded a scalar time-series xfc), 
i=l,...,N, at constant sampling time intervals Z\f=f, +r f, . We build m dimensional vectors in reconstructed state 
space as follows 



1 Published In: Zelinka I. (ed.): Perspective in Modern Prediction Methods, 44-48, Technical University of Brno, 
1998, ISBN: 80-214-1222-4 
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X(t x ) = {xfc ), x(t x + t), . . ., + (jw - 7)t)}, 



Ifc ) = {*(*,. ), xfc + t), . . ., x (t t +(m- 



(l) 



^ta ) = H f M \ A f M +*\~ -A*** +( m - J ^)l 



where M=N-(m-l)T, m is so-called embedding dimension and Tis an appropriate time delay (integer multiple 
of the sampling time At). Embedding theorem tells us that if the dynamical system evolves on a set of dimension 
D then sufficient condition for choosing embedding dimension m is m>2D. In experimental practice, D is 
unknown. Therefore, procedure of finding the embedding dimension is to increase m step by step and to compute 
a fractal dimension (or another invariant) for every value of m until the fractal dimension remains almost 
constant. At this moment we stop adding other components to the time delay vector. The time delay z can be 
chosen almost arbitrarily for an infinite amount of noise-free data. If the data are noisy and limited in number the 
value of T should guarantee an independence of the co-ordinates of the reconstructed state space. In this study, 
we select the time delay as the lag at which the autocorrelation function A(k) of the time series x(t : ) 



A{k) = 



where 



k=0,l,...,N-l 



(2) 



^ [x{t)-x\x{t 1+k -x)\ _ ;f 



attains the value of 1/e. We obtained T=63 days and T=4 days for temperatures and pressure time-series, 
respectively. The number of data points N=10957 for every of the sixteen time-series. For more sophisticated 
treatment on embedding procedure see, e.g., Sauer et al. (1991), Abarbanel et al. (1993). 

3. Correlation dimension 

We select correlation dimension Vas a first quantifier describing properties of a possible climate or 
weather attractor. Correlation dimension provides the number of active degree of freedom of the system. This 
dimension was introduced by Grassberger and Procaccia (1983a) because it is easily computable from 
experimental data. Correlation dimension is defined as follows 



r lnC m (r) 
V = lim lim — , 

r— >0 M— Jjj y 



(3) 



where C"(r) is the correlation integral 



(4) 



and 0(u) is Heavisede step function 
\0 if u < 



e{u)= 



1 if u > 
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1 1... 1 1 denotes an appropriate norm in the state space. The cut-off parameter T is introduced into denominator of 
Eq. (4) in order to eliminate an anomalous shoulder of correlation integral (see also Theiler 1986). In this 
contribution, we choose Tq=X. 

Extraction V directly from Eq. (3) is ineffective, therefore numerous estimators have been developed (for 
excellent review see Theiler 1990). We use the method suggested by Raidl (1996). This approach is based on 
searching for a "meaningful scaling region" in a plot of In C versus In r. Then V is estimated as a slope of the 
meaningful scaling region. 

Figure 1 shows the In C - In r diagram for the temperature time-series of Cheb station. Two scaling regions, 
denoted (A) and (B), are visible in this diagram. The slope of the region (A) oscillates between 20 and 30, 
whereas the slope of the region (B) converges quickly to the value of 1.7 (see Fig. 2). Analogous results were 
obtained for the other temperature time-series, with the exception of the time-series of the mountain stations 
(Praded, Lysa hora) - in these cases the slope of the region (B) converges to the value between 1 .9 and 2.0. The 
regions (B) disappear after subtraction an annual course of temperature from each daily value. 



Fig.l 
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Figure 1: In C - In r diagram for different embedding dimensions m (m=2, 4, 60) of the temperature 
time-series of Cheb station. 



Figure 2: Estimation of correlation dimension v as a function of embedding dimension m of the temperature 
time-series of Cheb station. 



Figure 3 shows typical correlation integrals of pressure time series. Only one scaling region is apparent in the 
InC-lnr diagram. Its slope oscillates between 20 and 30 (Fig. 4). 



Fig.3 




Figure 3: In C - In r diagram for different embedding dimensions m (m=2, 4, 60) of the pressure time-series 
of Cheb station. 

Figure 4: Estimation of correlation dimension v as a function of embedding dimension m of the pressure 
time-series of Cheb station. 
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4. Entropy 

An non-integer value of V does not guarantee that the underlying dynamics is chaotic. Therefore we 
proceed to estimation of Kolmogorov AT-entropy. It is well-known (Grassberger and Procaccia 1983b) that: 

• K=0 in ordered system, 



• K—> °° in purely random system, 



• < K < °° in chaotic (deterministic) system. 

We evaluate AT-entropy through A^-entropy which possesses the following properties: <K 2 , K 2 <K, and has the 
three above properties of K. K 2 is easy to compute by means of correlation integrals (Grassberger and Procaccia 
1983b): 



K 2 = limlim £"(/■), (5) 

m— >oo r— >0 

where 

K™=—\n C "ff, , (6) 
2 kAt C m+k {r) 



and k is sufficiently small integer number. In order to reduce fluctuations and to improve the statistics, we use 
averaging of Eq. (6) over five different values of k and approximate dependence of K 2 m on m by means of 
least-squares fit by the function 



K™(r) = K 2 +^-Y -In 
4LAtj-fl 



'm + 2l 



m 



L=5. 



(7) 



It stands to reason that A" 2 -entropy has to be evaluated in the same scaling region as v. Figure 5 depicts K 2 m as a 
function of m for temperature time-series of Cheb station. The other data provide similar results. The value of the 
error-doubling time T 2 is more illustrative of the predictability of the atmosphere than the value of A^-entropy. 
The error-doubling time is defined as follows 



Estimated values of T 2 are summarized in Tab. 2. Only the results of extrapolation from scaling regions (A) are 
mentioned because evaluations from scaling regions (B) give extremely high values of T 2 (for more information 
see Raidl 1995). 





Cheb 


Holesov 


Hradec 
Kralove 


Kuchafovic 


Mosnov 


Ruzyne 


Lysa hora 


Praded 


Temperature 


2.9±0.3 


1.8±0.2 


1.8±0.1 


2.1±0.2 


2.1±0.2 


2.1±0.3 


1.9±0.2 


2.1±0.1 


Pressure 


1.3±0.1 


1.1±0.1 


1.2±0.1 


1.2±0.1 


1.2±0.1 


1.1±0.1 


1.1±0.1 


1.1±0.1 



Table 1 : Error-doubling time T 2 [day] of the surface temperature and surface pressure time-series. 
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Figure 5: Plot of K 2 m as a function of embedding dimension m for the scaling region (A) of the temperature 
time-series of Cheb station. 



The ideas of non-linear time series analysis have been used to evaluate correlation dimension and 
A^-entropy of the atmospheric system. Two scaling regions have been revealed in the correlation integrals of the 
temperature time series. The first one possessed the gradual slope (V between 1.7 and 2.0) and the second one 
possessed rapid slope (v between 20 and 30). The first region is linked to the seasonal variability of weather 
because it disappeared from the In C - In r diagram after removing annual course of temperature from the 
time-series. Only one scaling region was detected in the correlation integral of the pressure time-series 
(v e<20, 30>). We conclude that the seasonal weather variability is governed by a low dimensional dynamical 
system, whereas the changes of weather from day to day are under the influence of high dimensional dynamics. 
The estimated values of K 2 greater than zero support the suggestion that the dynamics of weather system is 
chaotic. For that reason the weather is predictable over no more than short time interval ahead. Error-doubling 
time of the high dimensional system is less than 3 days. 
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